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ABSTRACT 


Independent Earth orientation measurements taken by the space-geodetic techniques of 
lunar and satellite laser ranging, very long baseline interferometry, and the Global Positioning 
System have been combined using a Kalman filter. The resulting combined Earth orientation 
series, SPACE2011, consists of values and uncertainties for Universal Time, polar motion, 
and their rates that span from September 28, 1976, to July 13, 2012, at daily intervals and is 
available in versions with epochs given at either midnight or noon. The space-geodetic 
measurements used to generate SPACE2011 have then been combined with optical astrometric 
measurements to form two additional combined Earth orientation series: (1) COMB2011, 
consisting of values and uncertainties for Universal Time, polar motion, and their rates that 
span from January 20, 1962, to July 13, 2012, at daily intervals and which are also available in 
versions with epochs given at either midnight or noon; and (2) POLE201 1, consisting of values 
and uncertainties for polar motion and its rate that span from January 20, 1900, to June 21, 2012, 
at 30.4375-day intervals. 
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INTRODUCTION 


Reference series of Earth orientation parameters (EOPs) obtained by combining 
independent measurements of the Earth’s orientation are generated annually at the Jet Propulsion 
Laboratory (JPL) in support of tracking and navigation of interplanetary spacecraft. This report 
describes the generation of the most recent such combined Earth orientation series: SPACE201 1, 
COMB201 1, and POLE201 1. Since the procedures used to generate these most recent series are 
similar to those used to generate previous such combinations, only a brief description of their 
generation is given here. Further details regarding the approach used at JPL to annually combine 
independent measurements of the Earth’s orientation can be found in Gross (1996, 2000) and 
Gross et al. (1998). 
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SPACE2011 


Data Sets Combined to Form SPACE2011 

SPACE2011 is a combination of independent space-geodetic measurements of the 
Earth’s orientation. Table 1 lists the space-geodetic series used in generating SPACE2011, 
giving their identifiers, the number of measurements from each series that were actually 
incorporated into SPACE201 1, and the time interval spanned by those measurements. Note that 
the University of Texas Center for Space Research (UTCSR) satellite laser ranging (SLR) 
Universal Time (UT) values were not used in generating SPACE2011 due to problems 
associated with separating this component of the Earth’s orientation from the effects of 
unmodeled forces acting on the satellite that cause the node of its orbit to drift (see Gross et al. 
1998, p. 217 for further discussion about this point). For similar reasons, the International Laser 
Ranging Service (ILRS) satellite laser ranging length-of-day (LOD) values have also not been 
used in generating SPACE201 1. 

Since it was desirable to combine only independent measurements of the Earth’s 
orientation, only those series listed in Table 1 were used, even though other space-geodetic 
series are available from other analysis centers. When more than one series determined by the 
same measurement technique was used, care was taken to make sure that the measurements 
themselves were not included more than once. In particular, polar motion measurements from the 
JPL Global Positioning System (GPS) series were only used until the start of the International 
Global Navigation Satellite Systems (GNSS) Service (IGS) combined series EOP(IGS) 95 P 01 
on January 1, 1995; polar motion measurements from the IGS combined series 

EOP(IGS) 95 P 0 1 were then used until the start of the IGS combined series EOP(IGS) 95 P 02 
on June 30, 1996; polar motion measurements from the IGS combined series EOP(IGS) 95 P 02 
were then used until the start of the accumulated IGS Solution Independent Exchange (SINEX) 
combined series EOP(IGS) 00 P 03 on February 26, 2000; and polar motion measurements from 
the accumulated IGS SINEX combined series EOP(IGS) 00 P 03 were used thereafter. Similarly, 
LOD measurements from the IGS combined series EOP(IGS) 95 P 02 were used until it ended 
on February 25, 2000 with LOD measurements from the accumulated IGS SINEX combined 
series EOP(IGS) 00 P 03 used thereafter. In addition, UT measurements from the National 
Oceanic and Atmospheric Administration (NOAA) International Radio Interferometric 
Surveying (IRIS) Intensive UT1 series were used until it ended on December 31, 1994; 
measurements from the United States Naval Observatory (USNO) National Earth Orientation 
Service (NEOS) Intensive UT1 series were then used until it ended on December 4, 2000; and 
measurements from the Goddard Spaceflight Center (GSFC) NEOS Intensive UT1 series were 
used thereafter. Finally, polar motion measurements from the UTCSR SLR series 
EOP(CSR) 96 L 01 were used until it ended on February 3, 1996 with measurements from the 
ILRS combined SLR series being used thereafter. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The Earth orientation series listed in Table 1 were first preprocessed by removing leap 
seconds from the Universal Time (UT1) values and, when necessary, by correcting the UT1 
values to be consistent with the extended definition of Greenwich Sidereal Time (GST) as 


2 



adopted by the International Earth Rotation and Reference Systems Service (IERS; IERS 1997, 
p. 149). Since most of the series listed in Table 1 were already consistent with the extended 
definition of GST, this correction needed to be applied to only the NOAA IRIS Intensive UT1 
series. 

Table 1 . Data Sets Combined to Fonn SPACE201 1* 


Data Set Name 

Data Type 

Analysis 

Center 

Reference 

Data Span 

Number 

LLR (JPL12M01; VOL, UTO) 
McDonald Cluster 

LLR 

JPL 

Williams et al. (2009) 

Oct. 5, 1976, to Jan. 5,2012 

747 

OCA 

LLR 

JPL 

Williams et al. (2009) 

April 7, 1984, to June 15, 2012 

1106 

Haleakala 

LLR 

JPL 

Williams et al. (2009) 

Feb. 10, 1985, to Aug. 11, 1990 

70 

Apache Point 

LLR 

JPL 

Williams et al. (2009) 

June 4, 2006, to July 2, 2012 

405 

UTCSR (CSR96L01) 
Lageos Polar Motion 

SLR 

UTCSR 

Eanes & Watkins (1996) 

Sep. 27, 1976, to Feb. 3, 1996 

2219 

ILRSA (30JUL12; Polar motion) 
ILRS Primary Combination 

SLR 

ASI-CGS 

Sciarretta et al. (2010) 

Feb. 4, 1996, to July 12, 2012 

5967 

DSN (JPL97R01; T, V) 
Califomia-Spain Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Nov. 26, 1979, to July 1 1, 2012 

1004 

California- Australia Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Oct. 28, 1978, to Sep. 30, 1997 

698 

NOAA (NOAA95R02) 
IRIS Intensive UT1 

VLBI 

NOAA 

Ray et al. (1995) 

April 2, 1984, to Dec. 31, 1994 

2394 

USNO (N9903) 

NEOS Intensive UT1 

VLBI 

USNO 

Eubanks et al. ( 1 999 ) 

Jan. 4, 1995, to Dec. 4, 2000 

1497 

NASA/GSFC (GSF2012a) 
NEOS Intensive UT1 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2012) 

Dec. 6,2000, to July 12, 2012 

3544 

NASA/GSFC (GSF2012a) 
Multibaseline 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2012) 

Aug. 4. 1979, to July 13,2012 

4434 

NASA/GSFC (GSFC 1 122) 
Westford-Fort Davis 

VLBI 

GSFC 

Gordon et al. (1999) 

June 25, 1981, to Jan. 1, 1984 

105 

Westford-Moj ave 

VLBI 

GSFC 

Gordon et al. ( 1 999) 

March 21, 1985, to Aug. 6, 1990 

18 

GPS (21APR04; Polar motion) 
Post-processed Flinn Analysis 

GPS 

JPL 

Heflin et al. (2004) 

June 10, 1992, to Dec. 31, 1994 

817 

GPS (IGS95P01; Polar motion) 
IGS Final Combined (ITRF93) 

GPS 

NRCan 

Kouba & Mireault (1997) 

Jan. 1, 1995, to June 29, 1996 

546 

GPS (IGS95P02; Polar motion) 
IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

June 30, 1996, to Feb. 25, 2000 

1336 

GPS (IGS00P03; Polar motion) 
IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to July 12, 2012 

4510 

GPS (IGS95P02; LOD) 

IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

Feb. 23, 1997, to Feb. 25, 2000 

1099 

GPS (IGS00P03; LOD) 

IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to July 12, 2012 

4520 


LLR, lunar laser ranging; JPL, Jet Propulsion Laboratory; VOL, variation of latitude; LIT, Universal Time; OCA, Observatoire de la Cote d’Azur; 
UTCSR, University of Texas Center for Space Research; SLR, satellite laser ranging; ILRS, International Laser Ranging Service; ASI, Agenzia Spaziale 
Italiana; CGS, Centro di Geodesia Spaziale; DSN, Deep Space Network; T, transverse; V, vertical; VLBI, very long baseline interferometry; NOAA, 
National Oceanic and Atmospheric Administration; IRIS, International Radio Interferometric Surveying; USNO, United States Naval Observatory; 
NEOS, National Earth Orientation Service; NASA, National Aeronautics and Space Administration; GSFC, Goddard Space Flight Center; GPS, Global 
Positioning System; IGS, International Global Navigation Satellite System (GNSS) Service; ITRF, International Terrestrial Reference Frame; NRCan, 
Natural Resources Canada; CODE, Center for Orbit Determination in Europe; IGN, Institut National de l'Information Geographique et Forestiere; 
SINEX, solution independent exchange; LOD, length of day. 


3 



Rotational variations caused by solid Earth and ocean tides were also removed from the 
UT1 values. The effect of the solid-Earth tides was removed by using the model of Yoder et al. 
(1981); the model of Kantha et al. (1998) was used to remove the effect upon UT1 of the ocean 
tides at the Mf, Mf', and Mm tidal frequencies. Since the Yoder et al. (1981) model already 
includes a contribution from the equilibrium ocean tides, only the Kantha et al. (1998) oceanic 
corrections to the Yoder et al. (1981) model were actually removed. Also note that the Kantha 
et al. (1998) model was used to remove the effect of ocean tides on only UT1, not on polar 
motion. Ocean-tide-induced polar motion variations were not removed from any of the polar 
motion observations. Finally, the only Earth orientation series listed in Table 1 that includes the 
effects of semidiurnal and diurnal ocean tides on the Earth’s orientation is the NOAA IRIS 
Intensive UT1 series. This series included these effects by adding to the released UT1 values the 
model of Herring (1993; also see Herring and Dong 1994). Hence, the same Herring (1993) 
model was used to remove them. 

On June 30, 1996, the IGS reference frame changed from the International Terrestrial 
Reference Frame ITRF93 to ITRF94; on March 1, 1998, it changed from ITRF94 to ITRF96; on 
August 1, 1999, it changed from ITRF96 to ITRF97; on February 27, 2000, it changed from 
ITRF97 to the IGS realization of ITRF97 known as IGS97; on December 2, 2001, it changed 
from IGS97 to the IGS realization of ITRF2000 known as IGS00; on January 11, 2004, it 
changed from IGS00 to a new IGS realization of ITRF2000 known as IGbOO; on November 5, 
2006 it changed from IGbOO to the IGS realization of ITRF2005 known as IGS05; on April 17, 
2011 it changed from IGS05 to the IGS realization of ITRF2008 known as IGS08 (Rebischung et 
al., 2011); and on October 7, 2012 IGS08 was replaced by the new IGS realization of ITRF2008 
kn own as IGb08 (Rebischung, 2012). These changes in reference frames potentially introduce 
discontinuities into the IGS combined Earth orientation series. However, the IGS Final 
combined series EOP(IGS) 95 P 02 and the accumulated IGS SINEX combined series 
EOP(IGS) 00 P 03 used here have had each of these segments aligned to the same IGS reference 
frame. Thus, to within the uncertainty in detennining the corrections required to align each 
segment, the IGS combined series EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 used here should be 
reasonably consistent with each other. They were, therefore, concatenated with one common set 
of bias-rate corrections being determined for them, as described below. 


Adjustments Made to Space-Geodetic Series Prior to Combination 

Prior to combining the series listed in Table 1 to form SPACE2011, series-specific 
corrections were applied for bias and rate, the stated uncertainties of the measurements were 
adjusted by multiplying them by series-specific scale factors, and outlying data points were 
deleted. Values for the bias-rate corrections and uncertainty scale factors were detennined by an 
iterative procedure wherein each series was compared, in turn, to a combination of all others. In 
order to minimize interpolation error (see Gross et al. 1998, pp. 223-225), the comparison of 
each series to its reference combination was done at the epochs of the measurements of that 
series by generating its reference combination using a Kalman filter that interpolates to (and 
prints its EOP estimates at) the exact epochs of those measurements. Also, both the bias-rate 
corrections and the uncertainty scale factors for all components of a given series were 
determined simultaneously in a multivariate approach using nonlinear weighted least squares. 
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Using a multivariate approach allows the correlations between the components to be taken 
into account when detennining the bias-rate corrections and uncertainty scale factors (see Gross 
et al. 1998, pp. 225). 

All the series listed in Table 1 were included in the iterative procedure. Details of the 
iterative procedure are described in Gross (1996, 2000) and Gross et al. (1998) and will not be 
repeated here. These details include: 

(1) Use of a reference series, SPACE2010 (Ratcliff and Gross, 2011), for initial bias-rate 
alignment; 

(2) Analysis of each data type in its natural reference frame; 

(3) Clustering of the McDonald LLR stations and, separately, the DSN very long baseline 
interferometry (VLBI) stations in California, Australia, and Spain; 

(4) Initial convergence on values for the series-specific bias-rate corrections and uncertainty 
scale factors prior to the start of outlier detection and deletion; and 

(5) Final convergence on these values after detecting and deleting all data outliers. 

At the end of the iterative procedure, relative bias-rate corrections have been detennined 
that make the series agree with each other in bias and rate; uncertainty scale factors have been 
determined that make the residual of each series have a reduced chi-square near one when 
differenced with a combination of all others; and outlying data points (those with residual 
values greater than four times their adjusted uncertainties) have been deleted. A total of 69 data 
points, or about 0.19% of those combined, were thus deleted. 

For GSFC NEOS Intensive and GSFC Multibaseline VLBI, the initial alignment to 
SPACE2010 was performed individually for each baseline containing at least 90 data points and 
spanning more than 1 year. Baselines containing less than one year of data or fewer than 90 data 
points were aligned as a group to SPACE2010. Subsequent steps of the iterative procedure were 
then conducted on the total intensive or multibaseline dataset as a whole. 

Starting April 17, 201 1 the IGS changed their procedure for estimating uncertainty in the 
GPS EOP values. This introduced a change in the magnitude of the polar motion X (PMX), polar 
motion Y (PMY), and LOD uncertainties of approximately a factor of 5. To account for this 
sudden change, bias, rate, and uncertainty scale factor values after April 17, 201 1 were estimated 
separately from those of earlier epochs. 

Finally, each series was placed within a particular IERS reference frame by applying 
to it an additional bias-rate correction that is common to all the series. This additional 
correction was determined by first combining all the series after applying to them the relative 
bias-rate corrections and uncertainty scale factors that had been determined for them as 
described above. This intermediate combination was then compared to the IERS combined 
Earth orientation series EOP(IERS) 08 C 04 (Bizouard and Gambis, 2011) during the interval 
1984-2011 in order to obtain the additional bias-rate correction required to make it, and 
therefore each individual series, agree in bias and rate with the IERS series. 

The total bias-rate corrections and uncertainty scale factors determined by the procedures 
outlined above are given in Table 2. The values of the bias-rate corrections given in Table 2 are 
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the sum of: (1) all the incremental corrections applied during the iterative procedure, (2) the 
corrections applied to initially align the series with each other, and (3) the additional, common 
correction applied in order to place each series within that particular IERS reference frame 
defined by the IERS combined Earth orientation series EOP(IERS) 08 C 04. The values of the 
uncertainty scale factors given in Table 2 are the products of all the incremental scale factors 
determined during the iterative procedure. The uncertainties of the bias-rate corrections given in 
Table 2 are the la standard errors of the incremental bias-rate corrections determined during the 
last iteration. There are no bias-rate entries in Table 2 for components that were either not used 
(e.g., the UTCSR SLR UT1 component) or not available (e.g., the NOAA IRIS Intensive polar 
motion components). Note that the same IERS rate correction is applied to all the data sets, 
including those such as the Westford-Mojave single baseline VLBI series, for which no relative 
rate correction could be determined. Therefore, the rate corrections given in Table 2 for those 
series for which no relative rate corrections could be determined are simply the IERS rate 
correction, but given in the natural reference frame of that series. In these cases, uncertainties 
for the rate corrections are not given. Also note that the entries for the bias-rate corrections in 
Table 2 for the IGS combined series EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 are the same 
prior to April 17, 2011. Since these entries were initially given within the same IGS reference 
frame, they were merged, and a common bias-rate correction was determined for them. 


Combined EOP Series: SPACE2011 

A Kalman filter was used to combine the series listed in Table 1 after the bias-rate 
corrections and uncertainty scale factors listed in Table 2 had been applied to them. The 
resulting combined Earth orientation series, SPACE2011, consists of values (Figure 1) and la 
standard errors (Figure 2) for polar motion, Universal Time, and their rates spanning September 
28, 1976, to July 13, 2012, at daily intervals; and it is available in versions for which the epochs 
are given at either midnight or noon. Leap seconds have been restored to UT1, and the effects 
of the solid Earth and ocean tides have been added back to the UT1 values using the same 
models for these effects that were originally used to remove them from the raw series: namely, 
the Yoder et al. (1981) model for the solid Earth tides and the Kantha et al. (1998) model for the 
ocean tides. However, semidiurnal and diurnal ocean tidal terms have not been added to and are 
therefore not included in the SPACE201 1 UT1 values. 

Figure 3 is a plot of the difference between the SPACE2011 polar motion, UT1, and 
LOD values and those of the IERS combined Earth orientation series EOP(IERS) 08 C 04. 
These two series are very consistent with each other, especially after January 1 , 2000 when the 
root-mean-square (rms) of their difference is only 0.054 milliarcseconds (mas) for the x- 
component of polar motion, 0.047 mas for the v-component, 0.020 milliseconds (ms) for UT1, 
and 0.018 ms for LOD. Prior to 2000, the difference between these two series exhibits greater 
variability and even some systematic behavior, particularly in the x-componcnt of polar motion. 
This systematic behavior is due to differences in the approaches used here and by the IERS to 
correct the bias and rate of the individual series before they are combined. 
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Table 2. Adjustments to Space-Geodetic Series ; 


Data Set Name 

Bias, 

mas 

Rate, mas/yr 

Uncertainty Scale Factor 

LLR (JPL12M01) 

VOL 

UT0 

VOL 

UT0 

VOL 

UT0 

McDonald Cluster 

-0.030 ±0.086 

0.150 ±0.088 

0.025 ±0.020 

-0.016 ±0.018 

1.474 

1.355 

OCA 

-0.041 ±0.041 

0.020 ± 0.033 

0.007 ± 0.006 

0.008 ± 0.004 

2.239 

1.897 

Haleakala 

—1.864 ± 1.178 

-1.542 ±0.743 

-0.342 ± 0.229 

-0.288 ±0.156 

2.744 

2.458 

Apache Point 

0.045 ±0.425 

-0.197 ±0.459 

-0.009 ± 0.025 

0.006 ± 0.027 

2.019 

2.056 

DSN (JPL97R01) 

T 

V 

T 

V 

T 

V 

Califomia-Spain Cluster 

-0.830 ±0.024 

-0.280 ± 0.059 

0.190 ±0.002 

0.172 ±0.005 

1.963 

1.417 

Califomia-Australia Cluster 

0.605 ±0.018 

0.476 ±0.053 

0.080 ± 0.008 

-0.003 ± 0.024 

1.445 

1.166 


NASA GSFC( 1122) T V 


Westford-Fort Davis 
Westford-Mojave 

8.764 ±3.786 
0.370 ± 1.054 

3.131 ±5.874 
-0.056 ± 1.647 

GSFC (GSF2012a) b 
G3GcKuRi 

PMX 

-0.009 ± 0.006 

PMY 

0.132 ±0.005 

UT1 

0.083 ±0.010 

HrRiWfWz 

0.572 ± 0.006 

-0.1 16 ±0.005 

0.085 ±0.010 

MoRiWfWz 

0.009 ± 0.006 

-0.283 ± 0.005 

-0.278 ±0.010 

Other networks 

0.004 ± 0.006 

0.009 ± 0.005 

0.017 ± 0.010 

NO A A (95R02) 
IRIS Intensive 

PMX 

PMY 

UT1 

0.195 ±0.023 

USNO (N9903) 
NEOS Intensive 

PMX 

PMY 

UT1 

-0.107 ±0.040 

GSFC (GSF2012a) b 
KkWz 

PMX 

PMY 

UT1 

-0.028 ± 0.025 

GnWz 



0.242 ± 0.025 

TsWz 



0.489 ± 0.025 

NyTsWz 



2.197 ±0.025 

KkSvWz 



0.621 ±0.025 

Other networks 



0.152 ±0.025 

UTCSR (96L01) 
Lageos 

PMX 

-0.394 ±0.010 

PMY 

0.730 ±0.009 

UT1 

ILRSA (30JUL12) 
Primary Comb. 

PMX 

0.025 ± 0.006 

PMY 

0.008 ± 0.006 

LOD 

GPS (21APR04) 

JPL Post-Flinn - 

PMX 

-0.270 ±0.016 

PMY 

0.137 ± 0.013 

LOD 

GPS (IGS95P01) 
Final Combined - 

PMX 

-0.136 ±0.056 

PMY 

0.191 ±0.054 

LOD 

GPS (IGS95P02) 
Final Combined 

PMX 

0.125 ±0.005 

PMY 

0.005 ± 0.005 

LOD 

-0.007 ±0.019 

GPS (IGS00P03) 
before 4/17/2011 
after 4/17/2011 - 

PMX 

0.125 ±0.005 
-0.012 ±0.288 

PMY 

0.005 ± 0.005 
0.785 ±0.278 

LOD 

-0.007 ±0.019 
-0.005 ± 1 .043 



T 

V 

T 


V 

0.859 ±0.381 0.284 ±0.588 

1.269 


0.909 

0.001 0.000 

2.364 


0.706 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

-0.037 ± 0.000 

0.057 ± 0.000 

-0.059 ±0.001 

1.569 

1.559 

1.699 

0.055 ± 0.000 

-0.1 16 ±0.000 

0.085 ±0.001 

1.569 

1.559 

1.699 

-0.072 ± 0.000 

-0.170 ±0.000 

-0.1 10 ±0.001 

1.569 

1.559 

1.699 

0.000 ± 0.000 

-0.001 ±0.000 

0.002 ±0.001 

1.569 

1.559 

1.699 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 



0.038 ± 0.007 



1.021 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 



0.041 ±0.007 



1.363 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 



0.003 ± 0.002 



1.653 



-0.003 ± 0.002 



1. 653 



-0.028 ± 0.002 



1. 653 



-0.135 ±0.002 



1. 653 



-0.037 ± 0.002 



1. 653 



-0.014 ± 0.002 



1. 653 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

0.101 ±0.004 

0.139 ±0.003 


0.742 

0.701 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

-0.000 ±0.001 

-0.007 ± 0.000 


3.065 

3.053 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

0.019 ±0.015 

-0.028 ±0.012 


1.312 

0.956 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

0.214 ±0.019 

0.321 ±0.019 


1.324 

0.998 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

-0.008 ± 0.000 

-0.001 ±0.000 

0.007 ± 0.002 

1.983 

1.583 

1.725 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

-0.008 ± 0.000 

-0.001 ±0.000 

-0.007 ± 0.002 

1.983 

1.583 

1.725 

0.002 ±0.015 

-0.042 ±0.015 

-0.520 ± 0.055 

12.617 

10.107 

8.518 


a ) Reference date for bias-rate adjustment is 1993.0. See Table 1 footnotes, mas, milliarcseconds; PMX, polar motion X; PMY, polar motion Y. 

b) GSFC VLBI Station IDs from ftp://cddis.gsfc.nasa.gov/vlbi/ivscontrol/ns-codes.txt 

G3 = NRAO 85-3; Gc = Gilmore Creek; Gn = NRAO 20m; Hr = HRAS at Fort Davis; Kk = Kokee Park; Ku = 9m at Kokee Park; Mo = 12m at 
Goldstone; Ny = Ny Alesund, Svalbard, Norway; Ri = Richmond, FL; Sv = Svetloe, Russia; Ts = Tsukuba 32m ; Wf = Westford; Wz = Wettzell. 
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A Combined Earth Orientation Series: SPACE201 1 
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Figure 1 . Plots of the x-component of polar motion (la), y-component of polar motion (lb), UT1-UTC (lc), and excess length-of-day (Id) as given by the 
combined Earth orientation series SPACE2011. The discontinuous changes in the plot of UT1-UTC are caused by the presence of leap seconds. Note that the 
UT1-UTC values displayed in (lc) include the tidal variations, whereas the excess length-of-day values shown in (Id) do not. 







Polar Motion X Error (mas) 


A Combined Earth Orientation Series: SPACE201 1 



Figure 2. Plots of the lo formal errors of the x-component of polar motion (2a), y-component of polar motion (2b), UT1-UTC (2c), and excess length-of-day 
(2d) as given by the combined Earth orientation series SPACE2011. The inserts within panels (2a), (2b), and (2c) show that component’s post-1984 
uncertainties on an expanded scale with the same units: milliarcseconds (mas) for polar motion, milliseconds (ms) for UT1-UTC. 
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Figure 3. Plots of the difference between the IERS combined Earth orientation series EOP(IERS) 08 C 04 and SPACE2011 formed by subtracting the SPACE2011 
values from those of the IERS series. The difference between the x-component of polar motion is shown in (3a), the difference between the y-component 
is shown in (3b), the difference between UT1-UTC is shown in (3c), and the difference between the excess length-of-day is shown in (3d). 




COMB2011 


COMB2011 extends SPACE2011 by additionally incorporating the optical astrometric 
polar motion and UT1 series that was detennined at the Bureau International de l’Heure (BIH) 
from an analysis of time and latitude observations by Li (1985; also see Li and Feissel 1986). 
This BIH optical astrometric series consists of values and uncertainties for polar motion and UT1 
that span from January 5.0, 1962 to December 31.0, 1981 at 5-day intervals. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The BIH optical astrometric series was first preprocessed by removing leap seconds 
from the UT1 values and by correcting the UT1 values to be consistent with the extended 
definition of GST, as adopted by the IERS (IERS 1997, p. 149). Rotational variations caused 
by solid Earth and ocean tides were also removed from the UT1 values. The same models that 
were used to remove the tidal effects from the series combined to form SPACE2011 were also 
used to remove them from the BIH series: namely, the Yoder et al. (1981) model for the solid 
Earth tides and the Kantha et al. (1998) model for the Mf, Mf, and Mm ocean tides. However, 
since the BIH UT1 measurements represent an average value over a 5 -day-long observation 
window, and since 5 days is a substantial fraction of the monthly and shorter-period tides, the 
amplitudes of these solid Earth and ocean tidal tenns were attenuated prior to their removal from 
the BIH UT1 measurements. (See Gross 1996, p. 8735 and Gross et al. 1998, pp. 226-227 for 
further discussion about this point.) 

Adjustments Made to BIH Series Prior to Combination 

The preprocessed BIH optical astrometric series was combined with the space-geodetic 
series that comprise SPACE201 1 after first: (1) correcting the BIH series to have the same bias, 
rate, annual terms, and semiannual terms as SPACE2011; (2) applying a constant multiplicative 
scale factor to the measurement uncertainties of the BIH series so that its residual, when 
differenced with SPACE2011, had a reduced chi-square of one; and (3) deleting those data 
points, if any, for which residual values were greater than four times their adjusted 
uncertainties. Due to software limitations associated with the need to correct the annual and 
semiannual terms of the BIH series, the above adjustments were determined separately for each 
component of the BIH series in a univariate approach rather than simultaneously in a 
multivariate approach as was done for the series combined to fonn SPACE201 1. The procedure 
used to determine these bias-rate and seasonal term corrections, uncertainty scale factors, and 
outlying data points has been described before (Gross 1996, pp. 8735-8738) and will not be 
repeated here. The annual and semiannual terms of the BIH series were adjusted in order to 
correct for systematic, seasonally varying effects that are known to be present in optical 
astrometric measurements. Since the values of both the BIH series and the SPACE2011 series 
are given at midnight, interpolation error (see Gross et al. 1998, pp. 223-225) is automatically 
minimized when differencing these two series for the purpose of determining the adjustments to 
be made to the BIH series. Tables 3 and 4 give the resulting uncertainty scale factors and values 
and la standard errors of the corrections to the bias, rate, annual, and semiannual terms thus 
determined for the BIH series. When determining these uncertainty scale factors and the 
corrections to the bias, rate, and seasonal terms, no outlying data points were detected. 
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Combined EOP Series: COMB2011 


A Kalman filter was used to combine the BIH series with the adjusted space-geodetic 
series that comprise SPACE2011 after first applying to the BIH series the uncertainty scale 
factors and corrections to the bias, rate, annual, and semiannual tenns that are given in Tables 3 
and 4. The resulting combined Earth orientation series, COMB201 1, consists of values (Figure 4) 
and la standard errors (Figure 5) for polar motion, Universal Time, and their rates that span from 
January 20, 1962, to July 13, 2012, at daily intervals and is available in versions with epochs 
given at either midnight or noon. Leap seconds have been restored to UT 1 , and the effects of the 
solid Earth and ocean tides have been added back to the UT 1 values using the same models for 
these effects that were originally used to remove them: namely, the Yoder et al. (1981) model for 
the solid Earth tides and the Kantha et al. (1998) model for the ocean tides. The full amplitude 
(i.e., no tidal terms attenuated) of the effects of the solid Earth and ocean tides at the epoch of the 
time tag were added back to the UT1 values. Semidiurnal and diurnal ocean tidal terms have not 
been added to and are, therefore, not included in the COMB201 1 UT1 values. 


Table 3. Adjustments to Bias, Rate, and Stated Uncertainty of Optical Astrometric Series* 


Data Set 

Bias, mas 



Rate, mas/yr 


Uncertainty Scale Factor 


PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

BIH 

-2.096 ±0.861 

-1.041 ±0.707 

-7.967 ± 1.081 

1.078 ±0.531 

1.130 ±0.433 

5.550 ±0.676 

1.830 

1.697 

1.950 

ILS 

-50.470 ± 2.384 

-2.906 ± 1.692 


0.444 ± 0.487 

-0.084 ±0.345 


2.184 

1.525 



*Reference date for bias-rate adjustment of BIH series is 1980.0. Reference date for bias-rate adjustment of ILS series is 1970.0. 


Table 4. Adjustment to Annual and Semiannual Terms of Optical Astrometric Series* 


Data Set 


Coefficient of Sine Term, mas 


Coefficient of Cosine Term, mas 


PMX PMY UT1 


PMX PMY UT1 


BIH 

annual -5.455 ± 1.061 -7.578 ± 0.892 5.305 ± 1.394 -2.693 ± 1.132 9.427 ± 0.933 -0.874 ±1.450 

semiannual 2. 189 ± 1.083 -0.481 ±0.907 -0.427 ± 1.421 0.982 ± 1.106 0.919 ±0.914 1.719 ±1.421 


ILS 

8.197 ±3.361 -10.512 ±2.381 

2.270 ±3.350 1.660 ±2.372 


annual -2.173 ±3.325 8.035 ± 2.352 

semiannual 0.178 ±3.334 8.813 ±2.359 


*Reference date for adjustment of annual and semiannual terms of BIH series is 1980.0. 
Reference date for adjustment of annual and semiannual terms of ILS series is 1970.0. 


12 



A Combined Earth Orientation Series: COMB201 1 






Figure 4. Plots of the x-component of polar motion (4a), y-component of polar motion (4b), UT1-UTC (4c), and excess length-of-day (4d) as given by the 
combined Earth orientation series, COMB2011. The discontinuous changes in the plot of UT1-UTC are caused by the presence of leap seconds. Prior to the 
introduction of leap seconds in 1972, Coordinated Universal Time (UTC) was adjusted by introducing step and rate changes designed to keep it close to UT1 
(e.g., IERS 1997, Table II-3), the effect of which is also readily apparent in (4c). Note that the UT1-UTC values displayed in (4c) include the tidal variations, 
whereas the excess length-of-day values shown in (4d) do not. 




















POLE2011 


No optical astrometric observations taken at the stations of the International Latitude 
Service (ILS) were used when creating the BIH optical astrometric series used in COMB2011 
(Li 1985; Li and Feissel 1986). The ILS polar motion measurements (Yumi and Yokoyama 
1980), which are based solely upon latitude observations made at the ILS stations are, therefore, 
independent of those comprising COMB2011 and have therefore been combined with them to 
form POLE201 1. Being based solely upon latitude observations, the ILS series contains no UT1 
measurements, but consists solely of polar motion measurements that span 1899.8-1979.0 at 
monthly intervals. Although no uncertainties are given with the individual polar motion values, 
the precision with which they have been detennined is estimated to be 10-20 mas (Yumi and 
Yokoyama 1980, p. 27). An initial uncertainty of 15 mas was, therefore, assigned to each of the 
ILS polar motion values. Since this assigned measurement uncertainty will be adjusted later, its 
initial value is arbitrary as long as it is not zero and serves merely as an a priori estimate to be 
used in the series adjustment procedure described below. 

The ILS series was combined with COMB201 1 to fonn POLE201 1 after: 

(1) correcting the ILS series to have the same bias, rate, annual, and semiannual 

terms as COMB201 1; 

(2) applying a constant multiplicative scale factor to the measurement uncertainties of 

the ILS series so that its residual, when differenced with COMB2011, had a 

reduced chi-square of one; and 

(3) deleting those data points, if any, whose residual values were greater than four 

times their adjusted uncertainties. 

Again, due to software limitations associated with the need to correct the annual and semiannual 
terms, these adjustments were determined separately for the x- and v-components of the ILS 
polar motion series by fitting a bias, rate, and these seasonal terms to the difference of the ILS 
series with COMB2011 during 1962.0 to 1979.0. The measurement uncertainties of the ILS 
polar motion values were adjusted by detennining and applying a scale factor that made the 
residual of this fit have a reduced chi-square of one. During this procedure to determine 
uncertainty scale factors and bias, rate, and seasonal term corrections, no outlying ILS data 
points were deleted since no data points had residual values greater than four times their adjusted 
uncertainties. Tables 3 and 4 (in the COMB2011 section) also give the resulting uncertainty 
scale factors and values and la standard errors of the corrections to the bias, rate, annual, and 
semiannual terms thus detennined for the ILS series. 

A Kalman filter was then used to combine the ILS series with the adjusted BIH and 
space-geodetic series that comprise COMB201 1, after applying to the ILS series the uncertainty 
scale factors and conections to the bias, rate, annual, and semiannual terms that are given in 
Tables 3 and 4. The resulting combined Earth orientation series, POLE201 1, consists of values 
(Figure 6a and 6b) and la standard errors (Figure 6c and 6d) for polar motion and its rate that 
span from January 20, 1900, to June 21, 2012, at 30.4375-day intervals. 
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Polar Motion X Error (mas) 


A Combined Earth Orientation Series: POLE201 1 
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Figure 6. Plots of the x-component of polar motion (6a), the ^-component of polar motion (6b), the 1g formal errors of the x-component of polar motion (6c), 
and the la formal errors of the j-component of polar motion (6d) as given by the combined polar motion series POLE201 1. 






DISCUSSION 


The Kalman filter that was used here to combine Earth-orientation measurements 
contains a stochastic model of the process that is used to propagate the state vector and its 
associated state covariance matrix forward in time to the epoch of the next measurement. For 
polar motion excitation, the stochastic model includes a random-walk term with equal noise 
forcing both the x- and v-components of excitation and originally having a white -noise power 
spectral density of 246.4 mas 2 /day (Morabito et al. 1988). This level of polar motion excitation 
process noise in the Kalman filter was increased to 739.2 mas 2 /day in order to better match the 
observed spectrum of polar motion excitation. Increasing the excitation process noise reduces 
the level of smoothing applied to the polar motion components of the propagated state vector and 
increases the covariance associated with those components. Thus, the SPACE2011, 
COMB201 1, and POLE201 1 polar motion and polar-motion excitation values are not as heavily 
smoothed as were those of SPACE2003, COMB2003, POLE2003, and earlier combinations 
produced at JPL; and the uncertainties assigned to the SPACE2011, COMB2011, and 
POLE201 1 polar motion and polar-motion excitation values are somewhat larger. 

The LOD components of the IGS Final combined EOP(IGS) 95 P 02 and EOP (IGS) 00 P 
03 series were incorporated into JPL’s combinations starting with SPACE2007 and COMB2007. 
In general, GPS LOD estimates are contaminated by orbital artifacts that must be removed prior 
to combination with other LOD or UT1 measurements (Chin et al., 2009). When determining 
the IGS Final combined LOD series, however, the IGS applies constraints that have the effect of 
removing these orbital artifacts. Therefore the IGS Final combined LOD series was incorporated 
into the Kalman filter as a true LOD measurement whose uncertainties are assumed to be 
“white” (i.e., random, uncorrelated in time, and with a flat power spectral density). Including the 
daily IGS Final combined LOD series in JPL’s combinations helps to compensate for less 
frequent UT1 measurements and was found to improve the agreement of JPL’s combined LOD 
estimates with independent atmospheric and oceanic angular momentum values, especially at the 
highest frequencies. 

Since a Kalman filter has been used in generating SPACE2011, COMB2011, and 
POLE20 11, the resulting polar motion and UT 1 values are smoothed to a degree that depends 
upon both the spacing between the measurements being combined (which determines how far the 
state vector and state covariance matrix must be propagated forward in time) and the 
uncertainties that have been assigned to the measurements. Since improvements to the 
observing systems, both in the hardware and software and in the number of systems, have led to 
increasingly precise determinations of the Earth’s orientation, and since the time resolution of 
the measurements has generally increased in concert with the measurement precision, the degree 
of smoothing applied to the SPACE2011, COMB2011, and POLE2011 values is a function of 
time, with the earlier values being more heavily smoothed than the more recent values. 

Daily EOP values are reported in SPACE2011 since the NOAA IRIS and GSFC NEOS 
Intensive UT1 values are given at daily intervals, as are the GPS and ILRSA combined SLR 
values (although gaps exist in each of these data sets). However, prior to the start of these data 
sets, the measurements combined to fonn SPACE2011 are given less frequently; therefore, the 
Kalman filter used to combine these measurements also interpolates them in order to produce a 
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series of equally spaced values. In order to be consistent with SPACE2011, daily EOP values 
are also reported in COMB2011 even though the BIH optical astrometric series used in 
COMB2011 is given at 5-day intervals. Thus, SPACE2011, COMB2011, and POLE2011 are 
equally spaced series of smoothed, interpolated Earth-orientation parameters. 

The combined Earth-orientation series SPACE2011, COMB2011, and POLE2011 are 
available from JPL’s Geodynamics and Space Geodesy Group via anonymous ftp: 
<ftp://keof.jpl.nasa.gov/keof/combinations/20 1 1> 
and upon request from the authors: Todd.Ratcliff@jpl.nasa.gov or Richard.Gross@jpl.nasa.gov. 
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Acronyms and Terms 


ASI Agenzia Spaziale Italiana 

BIH Bureau International de l’Heure 

CGS Centro di Geodesia Spaziale 

CODE Center for Orbit Determination in Europe 

COMB extends the SPACE series by additionally incorporating BIH optical astrometric 

measurements of polar motion and UT1 

DSN Deep Space Network 

EOP Earth orientation parameters 

GNSS International Global Navigation Satellite System 

GPS Global Positioning System 

GSFC Goddard Space Flight Center 

GST Greenwich Sidereal Time 

IERS International Earth Rotation and Reference Systems Service 

IGN Institut National de lTnformation Geographique et Forestiere 

IGS International Global Navigation Satellite Systems (GNSS) Service 

ILRS International Laser Ranging Service 

ILS International Latitude Service 

IRIS International Radio Interferometric Surveying 

ITRF International Terrestrial Reference Frame 

JPL Jet Propulsion Laboratory 

LLR lunar laser ranging 

LOD length-of-day 

NASA National Aeronautics and Space Administration 

NEOS National Earth Orientation Service 

NOAA National Oceanic and Atmospheric Administration 

NRCan Natural Resources Canada 

OCA Observatoire de la Cote d’Azur 

PMX polar motion X 

PMY polar motion Y 

POLE extends the COMB series by additionally incorporating ILS optical astrometric 

measurements of polar motion 

nns root mean square 
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SINEX 

SLR 

SPACE 


T 

USNO 

UT 

UT1 

UTC 

UTCSR 

V 

VLB I 
VOL 


Solution Independent Exchange 
satellite laser ranging 

a combination of independent space-geodetic measurements of the Earth's 
orientation 

transverse component of Earth orientation from single-baseline VLBI 

United States Naval Observatory 
Universal Time 

the principle fonn of Universal Time; conceptually, it is mean solar time 

Coordinated Universal Time 

University of Texas Center for Space Research 

vertical component of Earth orientation from single-baseline VLBI 
very long baseline interferometry 
variation of latitude 
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